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In this paper we give a complete analysis of the phase transitions in the mean-field Blume-Emery- 
Griffiths lattice-spin model with respect to the canonical ensemble, showing both a second-order, 
continuous phase transition and a first-order, discontinuous phase transition for appropriate values 
of the thermodynamic parameters that define the model. These phase transitions are analyzed both 
in terms of the empirical measure and the spin per site by studying bifurcation phenomena of the 
corresponding sets of canonical equilibrium macrostates, which are defined via large deviation prin- 
ciples. Analogous phase transitions with respect to the microcanonical ensemble are also studied via 
a combination of rigorous analysis and numerical calculations. Finally, probabilistic limit theorems 
for appropriately scaled values of the total spin are proved with respect to the canonical ensemble. 
These limit theorems include both central-limit-type theorems when the thermodynamic parameters 
are not equal to critical values and non-central-limit-type theorems when these parameters equal 
critical values. 
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I. INTRODUCTION 

The Blume-Emery-Griffiths (BEG) model [3] is an important lattice-spin model in statistical mechanics. 
It is one of the few and certainly one of the simplest models known to exhibit, in the mean-field approxi- 
mation, both a continuous, second-order phase transition and a discontinuous, first-order phase transition. 
Because of this property, the model has been studied extensively as a model of many diverse systems in- 
cluding He^-He^ mixtures — the system for which Blume, Emery, and Griffiths first devised their model 
[3] — as well as solid-liquid-gas systems [15, 21, 22], microemulsions [20], semiconductor alloys [16], and 
electronic conduction models [14]. On a more theoretical level, the BEG model has also played an impor- 
tant role in the development of the renormalization-group theory of phase transitions of the Potts model; see 
[13, 17] for details and references. 

As a long-range model with a simple description but a relatively complicated phase transition structure, 
the BEG model continues to be of interest in modern statistical mechanical studies. Our motivation for 
revisiting this model was initiated by a recent observation in [1, 2] that the mean-field version of the BEG 
model has nonequivalent microcanonical and canonical ensembles, in the sense that it exhibits microcanon- 
ical equilibrium properties having no equivalent within the canonical ensemble. This observation has been 
verified in [12] by numerical calculations both at the thermodynamic level, as in [1, 2], and at the level of 
equilibrium macrostates. In response to these earlier works, in this paper we address the phase transition 
behavior of the model by giving separate analyses of the structure of the sets of equilibrium macrostates 
for each of the two ensembles. Not only are our results consistent with the findings in [1, 2, 12], but also 
we rigorously prove for the first time a number of results that significantly generalize those found in these 
papers, where they were derived nonrigorously. For the canonical ensemble, full proofs of the structure of 
the set of equilibrium macrostates are provided. For the microcanonical ensemble, full proofs could not be 
attained. However, using numerical methods and following an analogous technique used in the canonical 
case, we also analyze the structure of the set of microcanonical equilibrium macrostates. 

The BEG model is a spin-1 model defined on the set {1, 2, n}. The spin at site j € {1, 2, n} is 
denoted by ujj, a quantity taking values in A = {—1, 0, 1}. The Hamiltonian for the BEG model is defined 
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In order to analyze the phase transition behavior of the model, we first introduce the sets of equihbrium 
macrostates for the canonical ensemble and the microcanonical ensemble. As we will see, the canonical 
equihbrium macrostates solve a two-dimensional, unconstrained minimization problem while the micro- 
canonical equilibrium macrostates solve a dual, one-dimensional, constrained minimization problem. The 
definitions of these sets follow from large deviation principles derived for general models in [7]. In the par- 
ticular case of the BEG model they are consequences of the fact that the BEG-Hamiltonian can be written 
as a function of the empirical measures of the spin random variables and that according to Sanov's Theorem 
the large deviation behavior of these empirical measures is governed by the relative entropy. 

We use two innovations to analyze the structure of the set of canonical equilibrium macrostates. The 
first is to reduce to a one-dimensional problem the two-dimensional minimization problem that characterizes 
these macrostates. This is carried out by absorbing the noninteracting component of the energy per particle 
function into the prior measure, which is a product measure on configuration space. This manipulation 
allows us to express the canonical ensemble in terms of the empirical means, or spin per site, of the spin 
random variables. Doing so reduces the analysis of BEG model to the analysis of a Curie-Weiss-type model 
[6] with single-site measures depending on (3. 

The analysis of the set of canonical equilibrium macrostates is further simplified by a second innovation. 
Because the thermodynamic parameter that defines the canonical ensemble is the inverse temperature (5, a 
phase transition with respect to this ensemble is defined by fixing the Hamiltonian-parameter K and varying 
/?. Our analysis of the set of canonical equilibrium macrostates is based on a much more efficient approach 
that fixes f3 and varies K. Proceeding in this way allows us to solve rigorously and in complete detail the 
reduced one-dimensional problem characterizing the equilibrium macrostates. We then extrapolate these 
results obtained by fixing (3 and varying K to physically relevant results that hold for fixed K and varying /3. 
These include a second-order, continuous phase transition and a first-order, discontinuous phase transition 
for different ranges of K. 

For the microcanonical ensemble, we use a technique employed in [1] that absorbs the constraint into 
the minimizing function. This step allows us to reduce the constrained minimization problem defining 
the microcanonical equilibrium macrostates to a one-dimensional, unconstrained minimization problem. 
Rigorous analysis of the reduced problem being limited, we rely mostly on numerical computations to 
complete our analysis of the set of equilibrium macrostates. Because the thermodynamic parameter defining 
the microcanonical ensemble is the energy per particle u, a phase transition with respect to this ensemble 
is defined by fixing K and varying u. By analogy with the canonical case, our numerical analysis of the 
set of microcanonical equilibrium macrostates is based on a much more efficient approach that fixes u and 
varies K. The analysis with respect to K rather than u allows us to solve in some detail the reduced one- 
dimensional problem characterizing the equilibrium macrostates. We then extrapolate these results obtained 
by fixing u and varying K to physically relevant results that hold for fixed K and varying u. As in the 
case of the canonical ensemble, these include a second-order, continuous phase transition and a first-order, 
discontinuous phase transition for different ranges of K. 

The contributions of this paper include a rigorous global analysis of the first-order phase transition in 
the canonical ensemble. Blume, Emery, and Griffiths did a local analysis of the spin per site to show that 
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their model exhibits a second-order phase transition for a range of values of K and that at a certain value of 
K a tricritical point appears [3]. This tricritical point has the property that for all smaller values of K, we 
are dealing with a first-order phase transition. Mathematically, the tricritical point marks the beginning of 
the failure of the local analysis; beyond this point one has to resort to a global analysis of the spin per site. 
While the first-order phase transition has been studied numerically by several authors, the present paper 
gives the first rigorous global analysis. 

Another contribution is that we analyze the phase transition for the canonical ensemble both in terms of 
the spin per site and the empirical measure. While all previous studies of the BEG model except for [12] 
focused only on the spin per site, the analysis in terms of the empirical measure is the natural context for 
understanding equivalence and nonequivalence of ensembles [12]. 

A main consequence of our analysis is that the tricritical point — the critical value of the Hamiltonian 
parameter K at which the model changes its phase transition behavior from second-order to first-order — 
differs in the two ensembles. Specifically, the tricritical point is smaller in the microcanonical ensemble 
than in the canonical ensemble. Therefore, there exists a range of values of K such that the BEG model 
with respect to the canonical ensemble exhibits a first-order phase transition while with respect to the mi- 
crocanonical ensemble the model exhibits a second-order phase transition. As we discuss in Section 5, 
these results are consistent with the observation, shown numerically in [12], that there exists a subset of the 
microcanonical equilibrium macrostates that are not realized canonically. This observation implies that the 
two ensembles are nonequivalent at the level of equilibrium macrostates. 

A final contribution of this paper is to present probabilistic limit theorems for appropriately scaled partial 
sums Sn = J2]=i with respect to the canonical ensemble. These limits, which follow from our work in 
Section 3 and known limit theorems for the Curie-Weiss model derived in [9, 10], include conditioned limit 
theorems when there are multiple phases. In most cases the limits involve the central-limit-type scaling n^/^ 
and convergence in distribution of Sn/n^^^ to a normal random variable. They also include the following 
two nonclassical cases, which hold for appropriate critical values of the parameters defining the canonical 
ensemble: 

Sn/n^^'^ — X, where P{X G dx} = const • exp[— const • x^] dx 

and 

Sn/n^^^ — ^ X, where P{X G dx} = const • exp[— const • x^] dx. 

As in the case of more compUcated models such as the Ising model, these nonclassical theorems signal the 
onset of a phase transition in the BEG model [6, Sect. V.8]. They are analogues of a result for the much 
simpler Curie-Weiss model [6, Thm. V.9.5]. 

The outline of the paper is as follows. In Section 2, following the general procedure described in [7], 
we define the canonical ensemble, the microcanonical ensemble, and the corresponding sets of equilib- 
rium macrostates. In Section 3, we outline our analysis of the structure of the set of canonical equilibrium 
macrostates. Because the proofs involve many technicalities, for ease of exposition we rely in this section 
on graphical arguments which, though not rigorous, convincingly motivate the truth of our assertions. The 
interested reader is referred to [18] for full details and complete rigor. In Section 4, we present new theoret- 
ical insights into, and numerical results concerning, the structure of the set of microcanonical equilibrium 
macrostates. In Section 5, we discuss the implications of the results in the two previous sections concern- 
ing the nature of the phase transitions in the BEG model, which in turn is related to the phenomenon of 
ensemble nonequivalence at the level of equiUbrium macrostates. Section 6 is devoted to probabiUstic Umit 
theorems for appropriately scaled sums Sji- 
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n. SETS OF EQUILIBRIUM MACROSTATES FOR THE TWO ENSEMBLES 

The canonical and microcanonical ensembles are defined in terms of probability measures on a sequence 
of probability spaces (A", The configuration spaces A" consist of microstates uj = (a;i, ...,a;n) with 
each ujj G A = {— 1,0,1}, and JF„ is the a-field consisting of all subsets of A". We also introduce the 
n-fold product measure P„ on VLn with identical one-dimensional marginals p = ^{6-i + 5o + 5i). 

In terms of the energy per particle hn,K defined in (1.1), for each n G IN, /? > 0, and K > the 
partition function is defined by 

Zni(3,K)= eyip[-nl3hn,K]dPn. 
J A" 

For sets B e J-'n,the canonical ensemble for the BEG model is the probability measure 

PnAK{B) = z;^^-^) • X exp[-n/3^„,i^] dPn. (2.1) 

For ueM, r>0, K>0 and sets B ^ Tn, the microcanonical ensemble is the conditional probability 
measure 

pn,r,K^^^ = I K,K €[u-r,u + r]} (2.2) 

^ Pn{B n {hn,K €[u-r,u + r]}} 
Pn{hn,K € [u - r,u + r]} 

As we point out after (2.4), for appropriate values of u and all sufficiently large n the denominator is positive 
and thus p^''^'^ is well defined. 

The key to our analysis of the BEG model is to express both the canonical and the microcanonical 
ensembles in terms of the empirical measure L„ defined for a; G A" by 

1 " 

ft - 

Ln takes values in "P = 'P(A), the set of probability measures on A. We rewrite hn^K as 

hn,K{oo) = — - K 

n \ n 



For iJ, define 



= J^y'^Lniuj,dy) - K (^J^yLn{uj,dy) 



(2.3) 



The range of this function is the closed interval [(1 — if) A 0, 1]. In terms of fx we express hn,K in the 
form 

We appeal to the theory of large deviations to define the sets of canonical equilibrium macrostates and 
microcanonical equilibrium macrostates. Since any fi E V has the form HiSi, where Hi > and 
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J2i=-i = 1» 7' can be identified with the set of probability vectors in M'^. We topologize V with the 
relative topology that this set inherits as a subset of M^. The relative entropy of /x G P with respect to p is 
defined by 

1 

i=-l 

Sanov's Theorem states that with respect to the product measures P„, the empirical measures L„ satisfy the 
large deviation principle (LDP) with rate function R{-\p) [6, Thm. Vin.2.1]. That is, for any closed subset 
FofVv/e have the large deviation upper bound 

limsup - logP„{L„ e F} < — inf R{p\p), 
n-»oo n neF 

and for any open subset G of P we have the large deviation lower bound 

limsup — log Pn{-^^n € G} > — inf R{p\p). 
n-*oo n neG 

From the LDP for the P„-distributions of we can derive the LDPs of L„ with respect to the two 

ensembles Pn,i3,K and P"'^'^. In order to state these LDPs, we introduce two basic thermodynamic func- 
tions, one associated with each ensemble. For /? > and > 0, the basic thermodynamic function for the 
canonical ensemble is the canonical free energy 

Vk{P) = - lim - log Zr,{(3,K). 

n— >oo 77, 

It follows from Theorem 2.4(a) in [7] that this limit exists for all /5 > and K > {) and is given by 

m{li)= inf{P(/x|p)+/3/K(/x)}. 

For the microcanonical ensemble, the basic thermodynamic function is the microcanonical entropy 

sk{u) = - inf {R{p\p) : /X G Mp) = u}. (2.4) 

Since R{p\p) > for all p, sk{u) £ [—oc, 0] for all u. We define dom s to be the set ofueM for which 
sk{u) > — oo. Clearly, dom sk coincides with the range of fx on V, which equals the closed interval 
[(1 — K) A0,1]. For u G dom sk and all sufficienriy large n the denominator in the second line of (2.2) is 
positive and thus the microcanonical ensemble P^'"''^ is well defined [7, Prop. 3.1]. 

The LDPs for L„ with respect to the two ensembles are given in the next theorem. They are conse- 
quences of Theorems 2.4 and 3.2 in [7]. 



Theorem 2.1. (a) With respect to the canonical ensemble Pn,i3,K' th^ empirical measures satisfy the 
LDP with rate function 

Ip,K{p) = R{p\p) + (ifK{p) - 'PKiP). (2.5) 

(b) With respect to the microcanonical ensemble P^'^'^, the empirical measures satisfy the LDP, in 
the double limit n ^ oo and r — > 0, with rate function 

P^i^{n) = [ -^(/^l^) + ^k{u) if f Kip) =U 

1 OO otherwise. 
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For n £ V and e > we denote by B{fi, e) the closed ball in V with center fi and radius e. If Ipdu.) > 0, 
then for all sufficiently small e > 0, inf^g^(^ IpifJ^) > 0. Hence, by the large deviation upper bound for 
Ln with respect to the canonical ensemble, for all € P satisfying Ip > 0, all sufficiently small £ > 0, 
and all sufficiently large n 

Pn,f3,K{Ln £ B{ii,e)} < exp -n(inf^^B{ix,6) Ipk^)) /2 

which converges to exponentially fast. Consequently, the most probable macrostates v solve I/3,k{i^) = 0. 
It is therefore natural to define the set of canonical equilibrium macrostates to be 

£/3,K = {lyeP: Ip,kM = 0} (2.7) 
= {v : u minimizes R{i^\p) + Pfxii^)}- 

Similarly, because of the large deviation upper bound for L„ with respect to the microcanonical ensemble, 
it is natural to define the set of microcanonical equilibrium macrostates to be 

5"'^ = {uGV: = 0} (2.8) 

= {v : v minimizes R{v\p) subject to fxiv) = u}. 

Each element v in Sp^K and £^'^ has the form v = z/_i(5_i + vq5q + vi5i and describes an equilibrium 
configuration of the model in the corresponding ensemble. For j = —1, 0, 1, Uj gives the asymptotic relative 
frequency of spins taking the value j. 

In the next section we begin our study of the sets of equiUbrium macrostates for the BEG model by 
analyzing Sp^x- 



m. STRUCTURE OF THE SET OF CANONICAL EQUILffiRIUM MACROSTATES 

In this section, we give a complete description of the set Sp^K of canonical equiUbrium macrostates for 
all values of (3 and K. In contrast to all other studies of the model, which fix K and vary (3, we analyze 

the structure of Ep^K by fixing (3 and varying K. As stated in Theorems 3.1 and 3.2, there exists a critical 
value of (3, denoted by (3c and equal to log 4, such that Sp^K has two different forms for [3 < /3c and for 
j3 > f3c. Specifically, for fixed (3 < f3c £p,K exhibits a continuous bifurcation as K passes through a critical 
value Kc^^ {(3), while for fixed (3 > £p,K exhibits a discontinuous bifurcation as K passes through a 

critical value Kc^^ {(3). In Section 5 we show how to extrapolate this information to information concerning 
the phase transition behavior of the canonical ensemble for varying (3: a continuous, second-order phase 
transition for all fixed, sufficiently large values of K and a discontinuous, first-order phase transition for all 
fixed, sufficiently small values of K. 

In terms of the uniform measure p = |(5-i + 6o + Si),v/e define 

ppidujj) = • exp(-/3a;|) p{dujj), (3.1) 

where Z{P) = exp(— p{dujj). The next two theorems give the form of £p^K for (3 < Pc and for 

P>/3c. 

Theorem 3.1. Define 13c = log 4 and fix P < Pc- Let pp be the measure defined in (3.1). The following 

results hold. 

(a) There exists a critical value Kc {(3) > such that 
(i)forK<KP{(3),£p^K = {pph 
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(ii) for K > Kc there exist probability measures v'^ = u'^{/3,K) and u~ = u~{(i,K) such 
that ^ ^ pp and = {^^j ^ 

(b) If we write as = u^-^S^i + Vq&q + Vi5i, then u~ = vfS^i + Uq&q + v^i&i. 

(c) and v~ are continuous functions of K > Kc and both pj^ and pp as 
Property (c) describes a continuous bifurcation in Sp^K ^.s K ^ {Kc^^ This analogue of a second- 

(2) 

order phase transition explains the superscript 2 on the critical value Kc The following theorem shows 
that for 13 > Pc the set Sp^x undergoes a discontinuous bifurcation as ii' ^ {Kc^\fi))^. This analogue of 
a first-order phase transtion explains the superscript 1 on the corresponding critical value Kc^\p). 

Theorem 3.2. Define fic = log 4 and fix P > Pc- Let pp be the measure defined in (3.1). The following 

results hold. 

(a) There exists a critical value Kc^^ {(5) > such that 

(i) forK <Ki^\p),8p,K = {pi}}; 

(ii) for K = Kc^\(3), there exist probability measures v'^ = v'^ ((3 , Kc^\(3)) and u~ = 
v~{(3, K^\f3)) such that ^ v~ ^ pp and £p,K = ■> ^ 

(Hi) For K > JTc^^ (/3), there exist probability measures v'^ = i''^((3, K) and v~ = y~((3, K) such 
that ^ v~ ^ pp and £p^K = {^^j ^ 

(b) If we write £ V as = v^id-i + Vq6o + Vidi, then v~ = I'ld-i + f(f(5o + I'lliSi. 

Since both v~^(l3,K) ^ pp and v^((3,K) / Pp for all K > K^\l3), the bifurcation in £p,K at 
K = Kc^\(3) is discontinuous. 

We prove Theorems 3.1 and 3.2 in several steps. In the first step, carried out in Section 3.1, we absorb 
the noninteracting component of the energy per particle into the product measure of the canonical ensemble. 
This reduces the model to a Curie- Weiss-type model, which can be analyzed in terms of the empirical means 

J2]=i i^j/n. The structure of the set of canonical equilibrium macrostates for this Curie-Weiss-type model 
is analyzed in Section 3.2 for (3 < fic and in Section 3.3 for f3 > I3c- Finally, in Section 3.4 we lift our 
results from the level of the empirical means up to the level of the empirical measures using the contraction 
principle, a main tool in the theory of large deviations. 



A. Reduction to the Curie- Weiss Model 



The first step in the proofs of Theorems 3.1 and 3.2 is to rewrite the canonical ensemble Pn,f3,K in the 
form of a Curie- Weiss-type model. We do this by absorbing the noninteracting component of the energy 
per particle hn,K into the product measure of Pn,i3,K- Defining Sn{uj) = Sj=i we write 



Pn,P,K(di^) = 



1 
1 

{z{0)T 

Zn{P,K) 



■ exp 



• exp 



• exp 




K 



n 



Pnidu) 



Y[ cxp{-Puj])p(duj] 
J j=i 

P^(d^). 
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In this formula Z{P) = ex.p(—PLOj)p{dujj) and is the product measure on A" with one-dimensional 
marginals pp defined in (3.1). 
We define 



/ exp 




/A" 





Since Pn,/3,k is a probability measure, it follows that 



and thus that 



Zn{P,K) 



■ exp 



Sniio) 



n 



(3.2) 



By expressing the canonical ensemble in terms of the empirical means Sn/n, we have reduced the 
BEG model to a Curie-Weiss-type model. Cramer's Theorem [6, Thm II.4.1] states that with respect to the 
product measure P^, Sn/n satisfies the LDP on M with rate function 



Jp{z) = SUp{tz - C/3{t)}. 

teiR 



In this formula cp is the cumulant generating function defined by 

cp{t) = log / exp{tLOi) pp{(kji) 

J A 



(3.3) 



(3.4) 



= log 



1 + 2e-l^ 



Jp is finite on the closed interval [—1, 1] and is differentiable on the open interval (—1, 1). This function is 
expressed in (3.3) as the Legendre-Fenchel transform of the finite, convex, differentiable function cp. By 
the theory of these transforms [19, Thm. 25.1], [6, Thm. VI.5.3(d)], for each z G (-1, 1) 



(3.5) 



From the LDP for 5„/n with respect to P^, Theorem 2.4 in [7] gives the LDP for Sn/n with respect to 
the canonical ensemble written in the form (3.2). 



Theorem 3.3. With respect to the canonical ensemble Pn,i3,K written in the form (3.2), the empirical means 
Sn/n satisfy the LDP on [—1, 1] with rate function 

If3,K = Mz) - fiKz'' - inf { J^(t) - ^Kt"]. (3.6) 

lEirt 

In Section 2 the canonical ensemble for the BEG model was expressed in terms of the empirical measures 
Ln- The corresponding set E^^k of canonical equihbrium macrostates was defined as the set of probability 
measures G P for which the rate function Ip^x in the associated LDP satisfies Ip^Ki^) = [see (2.7)]. 
By contrast, in (3.2) the canonical ensemble is expressed in terms of the empirical means Sn/n. We now 
consider the set S^^k of canonical equilibrium macrostates for the BEG model expressed in terms of the 
empirical means. The last theorem makes it natural to define as the set of z G [—1,1] for which the 
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rate function in that theorem satisfies Ip^xiz) = 0. Since z is a zero of this rate function if and only if z 
minimizes Jp{z) — fiKz"^, we have 

^p,K = {z & [—1, 1] : z minimizes Jp{z) — f3Kz^}. (3.7) 

As we will see in Theorem 3.8, each z G Sj^^x equals the mean of a corresponding measure in z/ G 
£l3,K- Thus, each z G Ep^K describes an equiUbrium configuration of the model in terms of the specific 
magnetization, or the asymptotic average spin per site. 

Although Ji3{z) can be computed expUcitly, the expression is messy. Instead, we use an alternative 
characterization of E^^k given in the next proposition to determine the points in that set. This proposition 
is a special case of a general result to be presented in [4]. 

Proposition 3.4. For z e M define 

Gp,K{z) = pKz^ - cp{2pKz). (3.8) 

Then for each /? > and K > 

min{J^(z) - pKz^} = uAn{Gf,,K{z)}. (3.9) 

In addition, the global minimum points of Jp{z) — j3Kz'^ coincide with the global minimum points ofGp^K- 
As a consequence, 

£f3^K = {z & M: z minimizes /^(z)}. (3.10) 
Proof. The finite, convex function f{z) = Cj3{2(5Kz) /2f3K has the Legendre-Fenchel transform 

We prove the proposition by showing the following three steps. 

1. sup,gjj{/(z) - z^/2] = sup|,|<i{zV2 - riz)}. 

2. Both suprema in Step 1 are attained, the first for some z & M and the second for some z E (—1, 1). 

3. The global maximum points of f{z)—z'^/2 coincide with the global maximum points of z^/2—f*{z). 
The proof uses three properties of Legendre-Fenchel transforms. 

1. For all zeM, f**{z) = {f*)*iz) equals f{z) [6, Thm. VI.5.3(e)]. 

2. If for some x e M and z e M we have z = f'{x), then f{x) + f*{z) = xz [19, Thm. 25.1], [6, 
Thm. VI.5.3(c)]. In particular, if z = x, then f{x) + f*{x) = x"^. 

3. If there exists x G (—1, 1) and y E M such that 

r(^) >r(x) + 2/(z-x)forall2G [-1,1], (3.11) 
then y = (/*)'(x) [19, Thm. 25.1]. Hence by properties 1 and 2 

f*{x) + f**{y) = f*{x) + f{y) = xy. 
In particular, if (3.11) is vaUd with y = x, then /(x) + f*{x) = x^. 
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Step 1 in the proof is a special case of Theorem C.l in [5]. For completeness, we present the straightfor- 
ward proof. Let M = sup^^julf (z) — z'^/2}. Since for any \z\ < 1 and x e M 

f*{z) + M>xz- f{x) +M>xz- x'^/2, 

we have 

f*{z) + M> sup{xz - x^/2} = z'^/2. 

It follows that M > z^/2 - f*{z) and thus that M > sup|^|<i{z^/2 - f*{z)}. To prove the reverse 
inequality, let N = sup|^|<]^{2;^/2 — f*{z)}. Then for any z & IR and jx| < 1 

z'^/2 + N>xz- x^/2 + N>xz- f*{x). 
Since f*{x) = oo for \x\ > 1, it follows from property 1 that 

z'^/2 + N> sup {xz - f*{x)} = f{z) 

\x\<l 

and thus that > sxi\)^^j[i{f{z) — z^ /2}. This completes the proof of step 1. 

Since f{z) ~ |z|as2;— >cxd,/(z) — attains its supremum over IR. Since /2 — /* (2;) is continuous 
and lim|2|^i(/*)'(z) = 00, z^/2 — f*{z) attains its supremum over [—1, 1] in the open interval (—1, 1). 
This completes the proof of step 2. 

We now prove that the global maximum points of the two functions coincide. Let x be any point in IR 
at which /(z) — 2^/2 attains its supremum. Then x = f'{x), and so by the second assertion in property 2 
/(•^) + f*{x) = x^- The point x lies in (—1, 1) because the range of f'{z) = </p{2(3Kz) equals (—1, 1). 
Step 1 now imphes that 

sup{/(z)-^V2} = f{x)-x''/2 

= xV2-r(x) = sup{zV2-r(^)}. 

\z\<l 

We conclude that z'^/2 — f*{z) attains its supremum at x € (—1, 1). 

Conversely, let x be any point in (—1, 1) at which 2:^/2 — f*{z) attains its supremum. Then for any 

-^€[-1,1] 

xV2-r(x)>zV2-r(z). 

It follows that for any 2; G [— 1, 1] 

r{z) > fix) + _ ^2)/2 > + _ a;). 

The second assertion in property 3 implies that f*{x) + f{x) = x^, and in conjunction with step 1 this in 
turn implies that 

sup{z^/2-riz)} = x^/2-rix) 

kl<l 

= f{x)-xy2 = sup{f{z)-zy2}. 
We conclude that f{z) — z'^/2 attains its supremum at x. This completes the proof of the proposition. ■ 
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Proposition 3.4 states that £p^K consists of the global minimum points of Gp^xiz) = PKz^—Ci3{2/3Kz). 
In order to simplify the minimization problem, we make the change of variables z — > z/2fiK in Gfj^x^ 
obtaining the new function 

FfsM^) = Gfs,K{z/2(3K) = cp{z). (3.12) 

Proposition 3.4 gives the alternative characterization of Sp^x to be 

£/3,K = e lR:w minimizes • (3.13) 

We use Fp^K to analyze £i3^k because the second term of Fp^x contains only the parameter /? while both 
terms in Gjs^x contain both parameters /? and K. In order to analyze the structure of S/^^k^ we take advantage 
of the simpler form of by fixing /3 and varying K. This innovation makes the analysis of Sgj^ much 
more efficient than in previous studies. Our goal is prove that the elements of Sp^K change continuously 
with K for all sufficiently small values of P [Thm. 3.1] and have a discontinuity at K = Kc^^ for all 
sufficiently large values of /3 [Thm. 3.2]. 

In order to determine the minimum points of Fp^x and thus the points in £/3,k< we study its derivative 

w 

Fp,KM = ^-Cp{w). (3.14) 

F'p j^{w) consists of a linear part w/2f3K and a nonlinear part c'p{w). According to Theorems 3.1 and 3.2, 
£l3^K exhibits a continuous bifurcation in K when /3 < /3c and a discontinuous bifurcation in K where 
j3c = log 4. As we will see in Sections 3.2 and 3.3, the basic mechanism underlying this change in the 
bifurcation behavior of £p^K is the change in the concavity behavior of (w) for P < Pc versus /3 > /3c, 
which is the subject of the next theorem. A related phenomenon was observed in [8, Thm. 1.2(b)] and in 
[11, Thm. 4] in the context of work on the Griffiths-Hurst-Sherman correlation inequality for models of 
ferromagnets; this inequality is used to show the concavity of the specific magnetization as a function of the 
external field. 

Theorem 3.5. Define f3c = log 4 and for (5 > j3c define 

WciP) = cosh"^ Qe^ - 4e-^^ > 0. (3.15) 

Then the following hold. 

(a) For (3 < Pc, c'piw) is strictly concave for w > 0. 

(b) For (3 > Pc> c^(^) strictly convex for < w < wdP) and c'p{w) is strictly concave for 
w > WciP). 

Proof, (a) We show that for all /3 < /3c, (^p{w) < for all w > 0. A short calculation yields 

[e-^(e^ _ e-^)][l - e-/^(e- + e^") - Se'^^] 

= [l+e-/3(e- + e-)]3 ^^'^^^ 

_ [2e-^ sinh w] [1 - 2e~l^ cosh w - Se'^/^] 
~ [1 + 2e-/^cosh-u;]3 

Since 2e~l^ sinhw and 1 -|- 2e~^ cosh are positive for w > c!'^{w) < for w > if and only if 

1 - 2e''^ cosh w; - 8e~^'' < for u; > 0. 
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The mequahty cosh u; > 1 for u; > implies that 

[1 - 2e-^coshu; - Se-^^] < [1 - 26"^ - Se'^^] = (1 - 4e-^){l + 2e-^) for all w > 0. 

Therefore, for all (3 < log 4, (^^'{w) < for tt; > 0. 

(b) Fixing /3 > /3c, we determine the critical value Wc{/3) such that c'p{w) is strictly convex for < ui < 
Wc{(3) and strictly concave for w > Wc{(3). From the expression for (^^{11]) in (3.16), cfp{w) > for > 
if and only if (1 — 2e~^ coshiu — Se"^'^) > for u; > 0. Therefore cfp{w) is strictly convex for 

< < cosh"^ Qe'^ - 4e~^^ . 

On the other hand, since c'^'(w) < for u; > if and only if (1 — 2e~^ coshty — 8e~^^) < for it; > 0, 
we conclude that (^p{w) is strictly concave for 

w > cosh~^ (^^^ ~ 4e~^^ . 
This completes the proof of part (b). ■ 

The concavity description of c/^ stated in Theorem 3.5 allows us to find the global minimum points of 

Ff^^K and thus the points in S^^k for all values of the parameters f3 and K. We carry this out in the next two 
sections, first for fi < j3c and then for (3 > [3c- In Section 3.4 we use this information to give the structure 
of the set £p^K of canonical equilibrium macrostates defined in (2.7). 



B. Description of £i3,k for /3 < /3c 

In Theorem 3.1 we stated the structure of the set Sj^^x of canonical equilibrium macrostates for the BEG 
model with respect to the empirical measures when P < Pc = log 4. The main theorem in this section, 
Theorem 3.6, does the same for the set Sfj^K, which has been shown to have the alternative characterization 

£/3,K = elR-.w minimizes Fp^Kiw)^ . (3.17) 

We recall that /^(u;) = iv^ /^j3K — Ci3{w), where cg is defined in (3.4). In Section 3.4 we use the fact 
that there exists a one-to-one correspondence between Ep^K and Sp^K to fully describe the latter set for all 
/3 < /3c and K > 0. 

According to Theorem 3.5, for P < /3c, c'p is strictly concave for u; > 0. As a result, the study of Sp^K is 
similar to the study of the equilibrium macrostates for the classical Curie-Weiss model as given in Section 
1V.4 of [6]. Following the discussion in that section, we use a graphical argument to motivate the continuous 
bifurcation exhibited by Sp^K for P < Pc- ^ detailed proof is given in Section 2.3.2 of [18]. 

Minimum points of Fp^x satisfy F'p ^i^) = 0, which can be rewritten as 

7/1 

= c'piw). (3.18) 



2PK 



Since the slope of the function w hh- w/2pK is 1/2PK, the nature of the solutions of (3.18) depends on 
whether 

4«»i^ or 0<^<4(0). 
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Define the critical point 

We use the same notation here as for the critical value in Theorem 3.1 because, as we will later prove, the 
continuous bifurcation in K exhibited by both sets E^^k and Ep^x occur at the same point k'^^ (/?) defined 
in (3.19). 

We illustrate the minimum points of -Fg,^- graphically in Figure 1 for /3 = 1. For three ranges of values 
of K this figure depicts the two components of F'^ the linear component w/2[3K and the nonlinear 

component c^(ii;). Figure 1(a) corresponds to < K < kP{(3). Since from (3.19) c^(0) = l/2pKP{/3), 

for < K < Kc (/3) the two components of Fp ^ intersect at only the origin, and thus Fp^^ has a unique 

(2) 

global minimum point at w = 0. Figure 1(b) corresponds to K = Kc {(i). In this case the two components 
of F'p ^ are tangent at the origin, and again Fp^x has a unique global minimum point at u; = 0. Figure 

1(c) corresponds to if > k'^^ (/?). For such K the global minimum points of Fp^K are symmetric nonzero 
points w = ±w{P, K), ■w{f3, K) > 0.1a addition, for K > Kc (P), w{(3, K) is a continuous function and 
as K ^ {K'i'\j3))^ , w{[3, K) converges to 0. As a result, we conclude that Sp^x exhibits a continuous 
bifurcation with respect to K. 



(a) (b) (c) 




Figures 1(a) and 1(c) give similar information as Figures IV. 3(b) and IV. 3(d) in [6], which depict the 
phase transition in the Curie- Weiss model. In these two sets of figures the functions being graphed are 
Legendre-Fenchel transforms of each other. 

Figures l(a)-(c) motivate the following theorem concerning the continuous bifurcation with respect to K 
exhibited by Sp^x in the BEG model. It is proved in Theorem 2.3.6 of [18]. The positive quantity K) in 
parts (b) and (c) of the next theorem equals w{f3, K)/2(5K, where K) is the positive global minimum 
point of Fp^K for K > K^ (/?). 

(2) 

Theorem 3.6. Define /3c = log 4 and fix (3 < (5c- Define £f3,K by (3.7) and the critical value Kc {(5) by 
(3.19). The following results hold. 

ia)ForK <KPif3),£p,K = {0}. 

(b) For K > Kc (/?), there exists a positive number z{f3^ K) such that £p^K = K)}, 

(2) 

(c) z{(3^K) is a strictly increasing continuous function of K > Kc (P), and z{P,K) as K ^ 

iKP{(3)r- 

This theorem completes our description of the continuous bifurcation exhibited by Sp^K for /3 < /3c- In 
the next section we describe the discontinuous bifurcation exhibited by Sp^K for P in the complementary 
region P > Pc- 
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C. Description of ^/3,if for (3 > j3c 

In Theorem 3.2 we gave the structure of the set Sp^K of canonical equilibrium macrostates for the BEG 
model with respect to the empirical measures when (3 > f3c. The main theorem in this section, Theorem 
3.7, does the same for the set S^^k, which has been shown to have the alternative characterization 

£q,k = 1 2^ ^IR-w minimizes F^,k(«^)| • (3.20) 

As in Section 3.2, Fg^^iw) = /A(3K — cp{w), where is defined in (3.4). In the next section we use 
the fact that there exists a one-to-one correspondence between Ep^x and Ep^x to fully describe the latter set 

for all (3 > (3,, and > 0. 

Minimum points of -F/3,k satisfy the equation 

F'p^xH = ^-c'f,iw)=0. (3.21) 

In contrast to the previous section, where for /3 < f3c c'^ is strictly concave for w > 0, part (b) of Theorem 
3.5 states that for (5 > (5c there exists Wc = Wc{(3) > such that is strictly convex for w G (0, Wc) 
and strictly concave for w > Wc- As a result, for (3 > Pc'WQ are no longer in the situation of the classical 
Curie-Weiss model for which the bifurcation with respect to K is continuous. Instead, for (3 > /3c, as K 
increases through the critical value Kc^\(3), Sp^x exhibits a discontinuous bifurcation. 

While the discontinuous bifurcation exhibited by £p^x for /3 > /3c is easily observed graphically, the full 
analytic proof is considerably more compUcated than in the case P < (3c- As in the previous section, we 
will motivate this discontinuous bifurcation via a graphical argument, referring the reader to Section 2.3.3 
of [18] for details. 

For (3 > (3cV/e divide the range of the positive parameter K into three intervals separated by the values 
Ki = Ki{(3) and K2 = K2{(3). Ki is defined to be the unique value of K such that the line w/2(3K is 
tangent to the curve at some positive point wi = wi{(3). The existence and uniqueness of Ki and wi are 
proved in Lemma 2.3.8 in [18]. K2 is defined to be the value of K such that the slopes of the line w/2(3K 
and the curve c|Jj at w = agree. Specifically, 

Figure 2 represents graphically the values of Ki and K2 for = 4. That figure exhibits Ki < K2; in 
Lemma 2.3.9 in [18] it is proved that this inequality holds for all (3 > (3c. 

In each of Figures 3-7, for fixed (3 > (3c and for different ranges of values of > 0, the first graph (a) 
depicts the two components of F'^ the Unear component w/2(3K and the nonlinear component c^. The 
second graph (b) shows the corresponding graph of F^^x- In these figures the following values of /3 were 
used: /3 = 4 in Figures 3, 5, 6, 7 and /? = 2.8 in Figure 4. 

As we see in Figure 3, for K G (0, Ki] the linear term intersects the nonlinear term at only the origin 
and thus Fp^x has a unique global minimum point at = 0. Since /? is fixed, the graph of the nonlinear 
term c'^^ also remains fixed. As the value of K increases, the slope of the Unear term w/2(3K decreases, 
leading to the discontinuous bifurcation in Ep^x with respect to K. 

The graph of Fp^x is depicted in Figure 4 for K G 00). We see that Fp^x has two global minimum 
points at u; = rk.w{(3, K), where w{(3, K) is positive. Therefore, by (3.20), for < K < Kx we have 
£0,K = {0} and for K>K2we have Sp^x = {±5(/3, K)}, where z{(3, K) = w{(3, K) /2(3K is positive. 

Now suppose that K G {Ki,K2). In this region, there exists w{(3,K) > such that Fp^x has three 
local minimum points at ui = and w = ±w(/3, K). As we see in Figure 5, for K slighty greater than 
Ki, Fp^xi^) < F/3,Kiw{P, K)); as a result, the unique global minimum point of Fp^x is w = 0. On the 
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1.5 




FIG. 2: Graphical representation of the values Ki and K2 for /? = 4. 




other hand, we see in Figure 6 that for K slightly less than K2, -f/3,i<:(0) > Fp^K{w{(3, K)); as a result, the 
global minimum points of Fj^^x are w = ±w{j3, K). As K increases over the interval (Ki,K2), F^ xi^) 
increases and Ffi^K{w{l3, K)) decreases continuously and consequently, as Figure 7 reveals, there exists 

a critical value K'i'\j3) such that F (1) (0) = F (1) Aw{[3, K)); as a result, the global minimum 

points of F ^(i).^x are li; = and w = ±w{/3, K). 

P^^c \P) 

In conclusion, we have the following picture: for Ki < K < kP{/3), S^^k = {0}; for K = KP{f3), 
£p,K = {0,±5(/3,if)}; and for K > kP{P), Sp^K = {±z{P,K)}. Lastly, since z{P,Ki^\p)) is 
positive, the bifurcation exhibited by Sp^x at X = Kc^^ (/?) is discontinuous. The same notation Kc^^ (/3) is 
used here to denote the critical value for the bifurcation exhibited by <?/3,a" as we used to denote the critical 
value for the bifurcation exhibited by Sp^K stated in Theorem 3.2. This is appropriate because the two 
critical values are equal. 

The discontinuous bifurcation exhibited by S/^^k for /? > /?c is described in the following theorem. This 

result corresponds to Theorem 2.3.7 in [18], where a complete proof is presented. The quantity z{l3,K) 
in parts (b)-(d) of the next theorem equals i/)(/3, K)/2[5K, where w{l3, K) is the positive global minimum 
point of Fp^K for K > K^'^^ (/3). 

Theorem 3.7. Define the set £/3,k i>y (3.7). For a fixed (3 > Pc = log 4, there exists a critical value k'^^ [P) 
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such that the the following hold. 

(a)ForK <kP{^),S0^k = {O}. 

(h)ForK = k''c\[3), £f3,K = {0,±z{p,K)} where 2{(3,K) > 0. 
{c)ForK> kP{I3), Sp^K = {±z{(3, K)} where zij3, K) > 0. 

(d) For K > Kc^\l3), z{(3,K) is a strictly increasing continuous function and z{P, Kc^\p)) > 0. 
Therefore, £p,K exhibits a discontinuous bifurcation. 

Together, Theorems 3.6 and 3.7 give a full description of the set Ss^k for all values of (3 and K. In the 
next section, we use the contraction principle to lift our results concerning the structure of the set Ei^^k up 
to the level of the empirical measures, making use of a one-to-one correspondence between the points in the 
two sets Ejj^K and Efj^x of canonical equihbrium macrostates. 

D. One-to-One Correspondence Between f /3,jf and Sj3,K 
We start by recalling the definitions of the sets £'^,x and Sjj^x- 

£p,K = {v eV -.v minimizes R{i'\p) + Pfxii^)} (3.23) 

and 

£p,K = {z z minimizes Jp{z) - fiKz^}. (3.24) 
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In the definition ofEp^K, R{iAp) the relative entropy of ji with respect to p = \{5-i+5q + 5\) and fxifJ-) 
is the function defined in (2.3). In the definition of £j3^K^ J/3 is the Cramer rate function defined in (3.3). 
We now state the one-to-one correspondence between the points in S^^k and the points in £i3,k- According 
to Theorems 3.6 and 3.7, S/^^k consists of either 1, 2 or 3 points. 

Theorem 3.8. Fix /? > and K > and suppose that £f3,K = {^A;}fc=i. r = 1,2 or 3. Define v^, k = 
1, r, to be measures in V with densities 

^{y) = eMtkV) ■ , (3.25) 



where is chosen such that y i^k{dy) = z^. Then for each k = 1, ...,r, tk exists and is unique, and £/3,k 
consists of the unique elements Vk, k = 1, . . . ,r. Furthermore, = 2PKzkfor k = 1, ...,r. 

The proof of the theorem depends on the following two lemmas. Both lemmas use the contraction 
principle [6, Thm. Vin.3.1], which states that for all z G [—1, 1] 

Jl3{z) = min : iieV, j^y ii{dy) = z| . (3.26) 

Lemma 3.9. ForP>OandK>0 

mm Ir{i^\pp) - f3K (^^y//(dy) M = min [jp{z) - f5Kz^] . 
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Proof. The contraction principle (3.26) implies that 

mml^R{li\pp)-pK(^J^yfi{dy)^ | 

= min ^min|i?(/i|pa) : ^ G V , j^y ij,{dy) = - ^Kz^ 

= mm{j,(.)-/3i^.^}. 
This completes the proof. ■ 

The second lemma shows that the mean of any measure v G Si^ k is an element of S/3,k- 
Lemma 3.10. Fix /3 > and K > 0. Given v € E-^^K' define z = Jj^y i^{dy). Then z G £i3,K- 
Proof. Since u G £p,K, is a global minimum point of — (3K{J^ y ^{dy)Y. Thus for all ;u G 7^ 

R{u\pp) -P^i^J^y u{dy)^ ' = R{u\pp) - (iK~z^ < R{p\pp) - PK (^J^y pi{dy)^ ' . 

In particular, this inequality holds for any jj, that satisfies ^^y ^{dy) = z. For such p,, the last display 
becomes 

Thus u satisfies 

The contraction principle (3.26) and Lemma 3.9 imply that 

Jf3{z)-PKz' = R{v\pp)-(5K(^j^yv{dy)^^ 

= ^^J^^R{p\Pt3) - (^J^yp{dy) 

= mm{Jg(z) - l3KzH. 
\z\<i '^^ ^ 

Therefore, z G £p,K, as claimed. This completes the proof. ■ 

We next prove Theorem 3.8. 
Proof of Theorem 3.8. A short calculation shows that for any /x G "P 

R{p\p) + PfKip) - ml{Ri,.\p)+PfKii^)} 

= R{p\p^) - (5K (^Jjpidy)^ - ind R{u\p^) - (5K (^j^V^dy) 
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Hence we obtain the following alternate characterization of Sp^K'- 

£/3,K = l^u -.u minimizes R{u\pfs) - f3K (^J y u{dy)j | . (3.27) 

For each Zk G £f),K, define the set 

^fe = |m e P : j^yix{dy) = Zfej . 

We first show for each A; = 1, .., r, is the unique global minimum point of R{fi\pi3) — f3K (/^ y ii{dy))'^ 
over A^. We then prove that 

i^l\^R{^i\p„) - PK Qji^idy)^^ = mi^l^R{fi\pfs) - (3K [Jj^dy) 

for all k,£ = 1, r. It will then follow that {i'k}k=i equals the set of global minimum points of R{iJL\pjj) — 
f3K (/^ y p,{dy))^ over the set A = Ufc=i ^fc- Finally, by showing that all the global minimum points of 
i?(/x|yO^) — f5K{J^ y ij,{dy))'^ lie in A, we will complete the proof that Sp^K = {i^k}k=i- If r = 2 or 3, then 
since y Vk{dy) = Zk, it is clear that if Zk ^ Z£, then Uk ^ v^. 

By Theorem Vlll.3.1 in [7], for each = 1, r, the point tk in the statement of Theorem 3.8 exists 
and is unique, 

Jp{zk) = R{uk\pp), (3.28) 

and R{iJi\pj3) attains its infimum over Ak at the unique measure Vk- Therefore, for each A; = 1, r, Vk is 
the unique global minimum point of R{p\pi3) — (3K (/^ y p{dy))^ over the set Ak- 
We next show that 

inf^l^R{p\pf,) - (3K (^J^ypidy)^^ = M^i^R{p\pp)-PK (^j^y pidy)^^ 

for all k,£ = 1, r. Since Zk, G £/3,K< Zk and Z(_ are global minimum points of Jp{z) — f3Kz'^. Thus by 
(3.28), we have 

inl\R{p\pp)-PK(j yp{dy)]\ = inf R{p\p^) - (5Kzl 

= Jfi{zk) - f3Kzl 
= Jp{z,)-(5Kzj 
= inf R{p\pf,)-pKz^, 

= m^lR{p\pp)-(5K(^Jjp{dy) 
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As a result, {i'k}k=i equals the set of global minimum points of R{p\pp) — (3K (/^ y p{dy))^ over the set 
A = \Jk=iAk. 

Lastly, we show R{p\pj3) — (3K{Jj^y p{dy))'^ attains its global minimum at points in A. Let cr be a 
global minimum point of R{p\pp) — (3K{J^y p{dy))^. By (3.27), this implies that a G Sp^K- Define 
C = IaV ^{dy). Then Lemma 3.10 implies that ( G £p^K and thus that C = for some k = 1, r. It 
follows that a e Ak C A for some A; = 1, r. 
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The last step is to prove that = 2(5Kzk for A; = 1, r. From definition (3.4), we have 

(^Sk) = I yn{dy) = Zk- 

J A 

In turn, the inverse relationship (3.5) implies that 

Therefore, since Zk G £i3,k, the definition (3.24) guarantees that Zk is a critical point of Jpiz) — f3Kz^. 
Thus, 

tk = 4{zk) = WKzk. (3.29) 
This completes the proof of Theorem 3.8. ■ 

E. Proofs of Theorems 3.1 and 3.2 

Theorem 3.1 gives the structure of the set Ep^K of canonical equilibrium macrostates, pointing out the 
continuous bifurcation exhibited by that set for f3 < (3c = log 4. The structure of Sp^x for (5 > /3c, given in 
Theorem 3.2, features a discontinuous bifurcation in K. The proofs of these theorems are immediate from 
Theorems 3.6 and 3.7, respectively, which give the structure of S/^^k for P < Pc and for /3 > /3c, and from 
Theorem 3.8, which states a one-to-one correspondence between Sp^x and Sp^K- 

Before proving Theorems 3.1 and 3.2, it is useful to express the measures pf^ and uj. in Theorem 3.8 
in the forms pp = pp-id-i + ppfi^o + pp,iSi and i^k = i^k,-iS-i + i^kfiSo + i^k,i^i, respectively. Since 
tk = 2f3Kzk, in terms of Zk £ S^^k we have 

e-P 1 e-P 



and 



Here 



In particular, Uk = P/3 when Zk = 0. 

We first indicate how Theorem 3.1 follows from Theorem 3.6. Fix /3 < Pc- The critical value K^\l3) 
in Theorem 3.1 coincides with the value K^\(i) in Theorem 3.6. For K < K^\l3), part (a) of Theorem 
3.6 indicates that S^^k = {0}; hence S/^^k = {pfs}- For K > Kc (P), part (b) of Theorem 3.6 indicates 
that Sf^^K = {^z{P, K)}, where z(P,K) > 0. It follows that the measures u'^ and v~ in part (a)(ii) of 
Theorem 3.1 are given by (3.25) with Zk = z{P, K) and Zk = —z{P, K), respectively. Since z{P, K) > 0, 
it follows that / / pp. Finally, part (c) of Theorem 3.6 allows us to conclude that u'^ and u~ are 
continuous functions of K > K^\p) and that both — > pp and ^ pp as K ^ {Kc^\p))'^. This 
completes the proof of Theorem 3.1. 

In a completely analogous way, Theorem 3.2, including the discontinuous bifurcation noted in the last 
line of the theorem, follows from Theorem 3.7. 



21 



In this section we have completely analyzed the structure of the set Sp^K of canonical equilibrium 
macrostates. In particular, we discovered that for f3 < (3c £0,K undergoes a continuous bifurcation at 
K = kP{P) [Thm. 3.1] and that for P > (3c £i3,K underg oes a discontinuous bifurcation at K = K^^ {(3) 
[Thm. 3.2]. We depict these bifurcations in Figure 8. While the second-order critical values K^\(3) are 
explicitly defined in Theorem 3.6, the first-order critical values k'^^ {(3) in the figure are computed numer- 
ically. The numerical procedure calculates k'^'^ {(3) for fixed values of (3 by determining the value of K 
for which the number of global minimum points of Gp^xiz) changes from one at z = to three at z = 
and z = ^z{(3, K), where z{(3, K) > 0. According to these numerical calculations for the discontinuous 

bifurcation, it appears that i^i^^ {(3) tends to 1 as /3 oo. However, we are unable to prove this hmit. 

In Section 5 we will see that Figure 8 is a phase diagram that describes the phase transitions in the 
canonical ensemble as /3 changes. We will also show that the nature of the bifurcations studied up to this 
point by varying K while keeping (3 fixed is the same if we vary (3 and keep K fixed instead. The latter 
situation corresponds to what is referred to physically as a phase transition; specifically, the continuous 
bifurcation corresponds to a second-order phase transition and the discontinuous bifurcation to a first-order 
phase transition. In order to substantiate this claim concerning the bifurcations and the phase transitions, 
we have to transfer our analysis of Sq^k from fixed (3 and varying if to an analysis of Sp^x for fixed K and 
varying (3. 

In the next section we study the BEG model with respect to the microcanonical ensemble. 
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FIG. 8: Bifurcation diagram for the BEG model with respect to the canonical ensemble 



rV. STRUCTURE OF THE SET OF MICROCANONICAL EQUILIBRIUM MACROSTATES 

In previous studies of the mean-field BEG model with respect to the microcanonical ensemble, results 
were obtained that either relied on a local analysis or used strictly numerical methods [1, 2, 12]. In this 
section we provide a global argument to support the existence of a continuous bifurcation exhibited by 
the set f of microcanonical equilibrium macrostates for fixed, sufficiently large values of u and for 
varying K. Specifically, for fixed, sufficiently large u E"^'^ exhibits a continuous bifurcation as K passes 
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through a critical value Kc^\u). The argument is similar to the one employed to analyze the canonical 
ensemble in Section 3. However, unlike the canonical case, where a rigorous analysis of the structure of 
the set Sp^K of canonical equilibrium macrostates was obtained for all values of /3 and K, the analysis of 
£U',K fQj. sufficiently large u and varying K relies on a mix of analysis and numerical methods. At the 
end of this section we summarize the numerical methods used to deduce the existence of a discontinuous 
bifurcation exhibited by 6'^'^ for fixed, sufficiently small u and varying K. In Section 5 we show how to 
extrapolate this information to information concerning the phase ttansition behavior of the microcanonical 
ensemble for varying u: a continuous, second-order phase transition for all sufficiently large values of K 
and a discontinuous, first-order phase transition for all sufficiently small values of K. 

We begin by recalling several definitions from Section 2. V denotes the set of probability measures with 
support A = {—1, 0, 1}; p denotes the measure + 6o + 5i) G for /x G "P 

1 



denotes the relative entropy of fi with respect to p; and fxifJ^) is defined by 

Mil) = jyix{dy)-Kiy^j^yix{dy)^^ 

= H-i) - K{p,i~ fi_if. 

For > we also defined the set of microcanonical equilibrium macrostates by 

f"'^ = {„eV: /"'^(z/) = 0} 

= {ly EV : u minimizes R{i^\p) subject to fxi^) = u}, 



(4.1) 



£u,K defined for > and u € doms^ = [(1 — K) AO, 1]. Throughout this section we fix 

u G domsft-- 

Determining the elements in f"'^ requires solving a constrained minimization problem, which is the 
dual of the unconstrained minimization problem associated with the set S^^k of canonical equilibrium 
macrostates defined in (2.7). In order to simplify the analysis of the set f we employ the technique used 
in [1] that reduces the constrained minimization problem defining S'^'^ to a one-dimensional, unconstrained 
minimization problem. For fixed K > and u G dom sk, we define 



1^u,K = {neV : fK{p) = u). 
For /Lt G Vu^K, let z = p\ — fi-i and q = pi + p-i. Since p G V^^k imphes that 

fxip) = [pi + P-\) - K{pi - p^if = u, 
we see that q = u + Kz^. Thus, for p eVuK, we have 



(4.2) 



i=-l 

"3 



2 

q + z 



log 



(q-z) 



log(gf + z) + 



+ (l-g)log[3(l-g)] + ii^log 
log(g - z) 



+ {l-q) log(l -q)-{q log 2 - log 3). 
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Setting q = u + Kz^, we define the quantity 

Ru,k{^) = ^ log(9 + z) + ^ logiq - z) (4.3) 
+ (!-(?) log(l -q)-iq log 2 - log 3) 

and the set 

Mu,K = {z e M : z = IJ.1 — for some /x G Vu^k}- (4.4) 

The derivation of Ru,k makes it clear that M.u,k C (—1, 1) is the domain of Ru,k- 
We next introduce the set 

^u,K ^ |~ g M.u,K '■ z minimizes Ru,k{z)}. 

The following theorem states a one-to-one correspondence between the elements of 8"^'^ and In [12], 
for particular values of u and K, numerical experiments show that ^"'^ consists of either 1, 2 or 3 points. 
Although we are not able to prove that this is valid for all u G dom sk and K > Q, because of our numerical 
computations we make it a hypothesis in the next theorem. 

Theorem 4.1. Fix K > and u G doms/^. Suppose ^"'^ = {zjfc}^^^, r = 1,..,3. Define Uk = 

J2i=-i ^k,iSi by the formulas 

. u + Kzl + Zk . u + Kzl — Zk 
^k,i- 2 ' ^k-i = 2 ' ^k,o = i- - t^Zk,l - J^Zk,-l- 

Then £^'^ consists of the distinct elements u^jk = 1, . . . ,r. 

Proof. Using the definition (4.2) of Vu^k, we can rewrite the set £^'^ of microcanonical equihbrium 
macrostates defined in (4.1) as 

gu,K _ g Vu^K : is a minimum point of R{fi\p)}. 

We show that for k = 1, .., r, fxii^k) = u and R{i'k\p) < R{l^\p) for all // G P^.k for which jj, ^ v^. 
From the definition of V]^ we have 

/if(ffc) = {i^k,! + J^k-i) - K{vk,i - T^k-if = {u + Kzl) - K^k = u. 
Therefore, Vk ^ T^u,K for all = 1, r. Since for all z^, zi G S'^'^ , k,£ = I, r 

R{vk\p) = Ru,K{zk) = Ru,K{ze) = R{Mp), 

it follows that R{i'f;\p) are equal for all k = 1, r. 

We now consider = Y^j^_i PiSi <E V^.k such that /v ^ for all A; = 1, r. Defining C, = pi — p-i, 
we claim that C / -Zfc for all A: = 1, r. Suppose otherwise; i.e. for some Zk 

m - ix_i = Q = Zk = i^k,! - ^k,-i- (4.5) 

But fi G Vu^K implies that fxip) = u = fxi^k) and thus that 

pi + /x_i = Vk^i + Uk-i. (4.6) 

Combining equations (4.5) and (4.6) yields the contradiction that = u^. Because C, ^ Zk for all k = 
1, ...,r, it follows that C, ^ £'^'^ and thus that Ru,K{zk) < Ru^k{0 for all A; = 1, ...,r. As a result, for 
k = 1, r we have 

R{vk\p) = RuM^k) < RuMO = R{p\p)- 



24 



We complete the proof by showing that if 7^ z^, then V}^ ^ v^. Indeed, if = v^, then for each 
choice of sign we would have Kz\ ± = Kz\ ± zi. Since this leads to the contradiction that = zi, the 
proof of the theorem is complete. ■ 

Theorem 4.1 allows us to analyze the set ^"'^ of microcanonical equilibrium macrostates by calculating 
the minimum points of the function Ru,k defined in (4.3). Define 

^u,k{z) = log{q + z) H — log(g - z) + {l-q) log(l - q), where q = u + Kz"^. 

With this notation (4.3) becomes 

RuA^) = ^uM^) - + ^2^) log 2 + log 3. 

This separation of Ru,k into the nonlinear component (pu,K and the quadratic component is similar to 
the method used in Sections 3.2 and 3.3 in determining the elements in the set £p,K- There we separated 
the minimizing function Fp^xiw) into a nonlinear component C/j(u;) and a quadratic component w'^/4(3K; 
minimum points of -F/3,a" satisfy F'^ ^{w) = c'g{2v) — w/2(3K = 0. Solving this equation was greatly 
facihtated by understanding the concavity and convexity properties of cp, which are proved in Theorem 3.5. 

Following the success of this method in studying the canonical ensemble, we apply a similar technique 
to determine the minimum points of Ru,k- We call a pair {u, K) admissible if n G dom sk- While an 
analytic proof could not be found, our numerical experiments show that there exists a curve K = C{u) 
in the {u, iC)-plane such that for all admissible {u, K) lying above the graph of this curve, ip'^ ^ is strictly 
convex on its positive domain. The graph of if = C{u) is depicted in Figure 9. We denote by the 
set of admissible (u, K) lying above this graph and by G_ the set of admissible (u, K) lying below this 
graph. Using a similar argument as in the proof of Theorem 3.6 for the canonical case, we conclude that 
for all {u, K) € the BEG model with respect to the microcanonical ensemble exhibits a continuous 
bifurcation in K; i.e., there exists a critical value Kc {u) > such that the following hold. 

1. ForK < iff^(u), ^"'-^ = {0}. 

2. For K > K^^ (u), there exists a positive number z{u, K) such that £^'^ = {±z(m, K)]. 

Combined with the one-to-one correspondence between the elements of ^"'^ and £'^'^ proved in Theo- 
rem 4.1, the structure of E"^'^ just given yields a continuous bifurcation in K exhibited by ^""-^ for {u, K) 
lying in the region G'^ above the graph of the curve K = C{u). Similar to the definition of the critical value 
Kc {(3) given in (3. 19) for the continuous bifurcation in K exhibited by £i3,k^ the critical value Kc {u) is 
the solution of the equation 

<x(0) = or <;^(0) = 2iflog2. 
Consequently, since (0) = l/u -I- 2ir[log(u/ (1 — u))], we define the second-order critical value to be 

if (2) (u) = = I (4 7) 

' ^ ^ 2 log 2 2.log (^)- ^ ^ ^ 

The derivation of this formula for Kc^^ {u) for the critical values of the continuous bifurcation in K exhibited 
by ^"'^ rests on the existence of the curve K = C{u), which in turn was derived numerically. However the 
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accuracy of (4.7) is supported by the fact that the graph of the curve Kc (n) fits the critical values derived 
numerically in Figures 2 and 3 of [12]. 

For values of {u, K) lying in the region G_ below the graph of the curve K = C{u), the strict convexity 
behavior of ip'^ ^ no longer holds. Therefore, numerical computations were used to determine the behavior 
of Ru,K for such (li, K), showing a discontinuous bifurcation in if in this region. Specifically, there exists 

a critical value K^'' (u) such that the following hold. 

1. For if < ifi^^(u), ^"'^ = {0}. 

2. For K = (u), there exists z{u, K) > such that ^"'^ = {0, ±z{u, K)}. 

3. For K > K^^^ (u), there exists z{u, K) >() such that 5"'^ = {±z{u, K)}. 

The critical values K^^^ (u) were computed numerically by determining the value of K for which the number 
of global minimum points of Ru,k{z) changes from one at 2 = to three at 2; = and z = ±z{u, K), 
z{u,K)>Q. 

The results of this section are summarized in the bifurcation diagram for the BEG model with respect 
to the microcanonical ensemble, which appears in Figure 9. In the next section we will see that Figure 9 is 
a phase diagram that describes the phase transition in the microcanonical ensemble as u changes. In order 
to substantiate this, we have to transfer our analysis of 6""'^ from fixed u and varying if to an analysis of 
fQj. fixed K and varying u. 
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FIG. 9: Bifurcation diagram for the BEG model with respect to the microcanonical ensemble 



V. COMPARISON OF PHASE DIAGRAMS FOR THE TWO ENSEMBLES 

We end our analysis of the canonical and microcanonical ensembles by explaining what our results 
imply concerning the nature of the phase transitions in the BEG model. These phase transitions are defined 
by varying (3 and u, the two parameters that define the ensembles. As we will see, the order of the phase 
transitions is a structural property of the phase diagram in the sense that it is the same whether we vary K or 
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(3 in the canonical ensemble and if or u in the microcanonical ensemble while keeping the other parameter 
fixed. 

Before doing this, we first review one of the main contributions of the preceding two sections, which is to 
analyze the bifurcation behavior of the sets E^^k and 5"'^ of equiUbrium macrostates with respect to both 
the canonical and microcanonical ensembles. Figure 8 summarizes the canonical analysis and Figure 9 the 
microcanonical analysis. The figures exhibit two different values of K called tried tical values and denoted 
by and K^'^^°. As we soon explain, at each of these values of K the corresponding ensemble 

changes its behavior from a continuous, second-order phase transition to a discontinuous, first-order phase 
transition. 

For the canonical ensemble, the tricritical value in Figure 8 is given by 

^canon ^ Kf^^^^) = (log 4) ^ 1.0820, 

(2) 

where Kc {(3) is defined in (3.19). With respect to the microcanonical ensemble, the tricritical value 

('2) 

j;^micro |g jj^g valuc of K at which the curves K = C{u) and Km {u) shown in Figure 9 intersect. From 
the numerical calculation of the curve K = C{u) we obtain the following approximation for the tricritical 
value Kg"'": 

if^'^'"" fti 1.0813. 

These values of K^°^ and ifJJ-'"" agree with the values derived in [1] via a local analysis and numerical 
computations. 

We first illustrate how our analysis of Sp^K in Theorems 3.1 and 3.2 for fixed (3 and varying K yields 

a continuous, second-order phase transition and a discontinuous, first-order phase transition with respect to 
the canonical ensemble. These phase transitions are defined for fixed K and varying j3, the thermodynamic 
parameter that defines the ensemble. In order to study the phase transition, we must therefore transform the 
analysis of S/^^k for fixed (3 and varying if to an analysis of the same set for fixed K and varying (3. After 
we consider the microcanonical phase transition in an analogous way, we will focus on the region 

k:^^'° 1.0813 <K < 1.0820 ^ K^°''. 

As we will point out, the fact that for K in this region the two ensembles exhibit different phase transition 
behavior — discontinuous for the canonical and continuous for the microcanonical — is closely related to 
the phenomenon of ensemble nonequivalence in the model. 

We begin with the continuous phase transition for the canonical ensemble. Figure 8 exhibits a mono- 

(2) 

tonically decreasing function K = Kc (J3) for < /? < /3c = log 4. Inverting this function yields a 
monotonically decreasing function (3 = (K) for K > Kf^f"™ = K^^\pc) ~ 1.0820. Consider for 
fixed K > K^"""" and small 6 > 0, values of /3 G {I3P{K) - 6, l3f\K) + 5). Our analysis of Ep^x in 
Theorem 3.1 shows the following. 

(2^ (2^ 

• For (3 G {13c {K) — 5, (3c {K)] the model exhibits a single phase pp. 

• For /3 G {pP (K), (K) + S) the model exhibits two distinct phases K) and W (/?, K). 

We claim that for fixed K > if^?"™ this is a second-order phase ttansition; i.e., as /3 — > {Pc^\k))~^ we 

have K) pp and K) pp. To see this, we recall from Figure 1(b) that for /? = P^\k) 
the graph of the hnear component wjipK of F'p p-{w) is tangent to the graph of the nonUnear component 
c'p{w) of F'p ^{w) at the origin. This figure was used in Section 3.1 to analyze the structure of the set Sp^K 
[Thm. 3.6]. Since both components of F'g j^{w) are continuous with respect to /3, a perturbation in /? yields 
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a continuous phase transition in Ep^x and thus in Sp^K- A similar argument shows that each of the double 

phases K) and K) are continuous functions of (3 for (3 > (3).' (K). 

We now analyze the discontinuous phase transition for the canonical ensemble in a similar way. Figure 
8 exhibits a monotonically decreasing function K = Kc^^ for P > f3c = log 4. Inverting this function 
yields a monotonically decreasing function (3 = /3c^\k) for K < 1.0820. Consider for fixed 

K < and small 5 > 0, values of /3 G {pi^^ {K) - 5, pi^^ (K) + (5). Our analysis of Sp,K in Theorem 

3.2 shows the following. 

• For P G {f3c^^ {K) — 6, /3c^-* (K)) the model exhibits a single phase pp. 

• For P = Pc^\k) the model exhibits three distinct phases pp, u~^{P, K), and v~{P, K). 

• For /3 G [pi^^ (K), pi^^ (K) + S) the model exhibits two distinct phases K) and u'ip, K). 

We claim that for fixed K < this is a first-order phase transition; i.e., as P —>■ {Pc^\k))~^, 

we have for each choice of sign v^{P, K) v^{Pc^\k), K) / pp. To see this, we recall from Figure 

7(a) that for P = P^\k) the graph of the linear component w/2PK of Fpj^{w) intersects the graph 
of the nonUnear component cfp{w) of Fg f^{w) in five places such that the signed area between the two 
graphs is 0. This results in three values of w that are global minimum points of Fp^K', namely, w = 
0,w{p, K), —w{P, K) [Thm. 3.7]. These three values of w give rise to three values of z = w/2PK 
that constitute the set £p^K for P = Pc {K). Since both components of F'^ ^{w) are continuous with 
respect to P, a perturbation in P yields a discontinuous phase transition in Sp^x and thus £p,K- A similar 
argument shows that each of the double phases iy'^{P, K) and K) are continuous functions of P for 
/3>#(if). 

The phase transitions for the microcanonical ensemble are defined for fixed K and varying u, the ther- 
modynamic parameter defining the ensemble. Therefore, in order to study these phase transitions, we must 
transform the analysis of f""'^ done in Section 4 for fixed u and varying K to an analysis of the same set 
for fixed K and varying u. This is carried out in a way that is similar to what we have just done for the 
canonical ensemble. In particular, we find that for K > ^ 1.0813 the BEG model with respect to 

the microcanonical ensemble exhibits a continuous, second-order phase transition and that for K < K^''^° 
the model exhibits a discontinuous, first-order phase transition. 

We now focus on values of K satisfying K^^^^° < K < K^^f°'^. As we have just seen, for such K 
the two ensembles exhibit different phase transition behavior: for i^"™"" < K the microcanonical ensem- 
ble undergoes a continuous, second-order phase transition while for K < K^°^ the canonical ensemble 
undergoes a discontinuous, first-order phase transition. This observation is consistent with a numerical cal- 
culation given in Figure 10 showing that for a fixed value of K G K^^^™) there exists a subset of 
the microcanonical equilibrium macrostates that are not realized canonically [12]. As a result, for this value 
of K the two ensembles are nonequivalent at the level of equilibrium macrostates. 

Figures 10(a) and 10(b) exhibit, for a range of values of n and P, the structure of the set 5"'^ of micro- 
canonical equiUbrium macrostates and the set Sp^x of canonical equiUbrium macrostates for K = 1.0817. 
This value of K hes in the interval {K^^""", Ktri"°") (1.0813, 1.0820). Each equilibrium macrostate in 
£u,K empirical measure having the form 

u = uiSi + uoSq + u^iS-i. 

In both figures the solid and dashed curves can be taken to represent the components ui and u-i. The 
components ui and in the microcanonical ensemble are functions of u [Fig. 10(a)] and in the canonical 
ensemble are functions of P [Fig. 10(b)]. Figures 10(a) and 10(b) were taken from [12]. 
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HG. 10: Structure of (a) the set f "'^ and (b) the set Efi^K for K = 1.0817. 



Comparing the two figures reveals that the ensembles are nonequivalent for this value of K. Specifically, 
because of the discontinuous, first-order phase transition in the canonical ensemble, there exists a subset of 
V that is not realized by S-fi^K for any /3 > 0. On the other hand, since the set E"^'^ of microcanonical 
equilibrium macrostates exhibits a continuous, second-order phase transition, the subset of V not realized 
canonically is realized microcanonically. As a result, there exists an inequivalence of ensembles at the 
level of equilibrium macrostates. The reader is referred to [12] for a more complete analysis of ensemble 
equivalence and nonequivalence for the BEG model. 



VI. LIMIT THEOREMS FOR THE TOTAL SPIN WITH RESPECT TO Pn,/3,K 

In Section 3.1, we rewrote the canonical ensemble Pn,ji,K in terms of the total spin Sn and thus reduced 
the analysis of the set E'^^k of canonical equilibrium macrostates to that of a Curie-Weiss-type model. 
We end this paper with Umit theorems for the P„ i^-distributions of appropriately scaled partial sums 
Sn = J2]=i '^j- represents the total spin in the model. Since Sn/n = J\y Ln{dy), the limit theorems 
for Sn are also limit theorems for the empirical measures L„. The new theorems follow from limit theorems 
for the Curie-Weiss model proved in [9] and [10]. 

Recall the function Gf^^K defined in (3.8) as 

Gp^K{z) = mz'-cp{2[5Kz), (6.1) 

where cp is defined in (3.4). Proposition 3.4 characterizes the set of canonical equilibrium macrostates £p^K 
by the formula 

^p.K = {z G IR : z minimizes Gp^xiz)}. 

In [9] and [10], it is proved that the limits in distribution of appropriately scaled partial sums 5„ in the 
Curie-Weiss model are determined by the minimum points of an analogue of Gp^K- As defined in [10, eqn. 
(1.6)], this analogue is 



Gp{z) = ^Pz^-cp{pz). 
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Carrying out a similar analysis of the minimum points of G^^k yields limit theorems for appropriately 
scaled partial sums Sn for the BEG model. The Umit theorems for the BEG model are proved exactly as in 
the Curie-Weiss case. 

The function Gp^x is real analytic. Hence for each global minimum point z = z{P, K) G £j3,K, there 
exists a positive integer r = r{z) such that > and 

Gp,K{z) = Gp,K{z) + + 0{{z - ~zr+^) asz^ ~z. 

We call r{z) the type of the minimum point z. This concept is well-defined since is real analytic and 
5 is a global minimum point. 

Because the limiting distributions for the scaled partial sums depend on the type of the minimum points 
z, we now classify each of the points in Sf^ x by type. This is done in Theorem 6.1 for (5 < f3c and K > 
in which case Ep^K exhibits a continuous bifurcation, and in Theorem 6.2 for (3 > /3c and i^T > 0, in which 
case £p^K exhibits a discontinuous bifurcation. The associated limit theorems are given in Theorems 6.3 and 
6.4. In all cases but one [Thm. 3.6(b)] the type of each of the minimum points is 1. When the type is 1, the 
associated limit theorems are central-Umit-type theorems with scalings n^/^. If = {0}, then Sn/n^^^ 
converges in distribution to an appropriate normal random variable, and if Sf^^x consists of multiple points, 
then Sn/n^^"^ satisfies a conditioned central-limit-type theorem. On the other hand, when K = K^iP), the 
type of the minimum point at is r = 2 or r = 3 depending on whether (3 < f5c or (3 = (3c- The associated 
limit theorems have non-central-limit scalings n^~^^^'^, and 

Pn,i3,K{Sn/n'^~^^^^ € dx} =^ const • exp[— const ■ x^^] dx. 

These non-classical limit theorems signal the onset of a phase transition [6, Sect. V.8]. 

We first consider l3 < (3c = log 4. According to Theorem 3.6, in this case there exists a critical value 

with the following properties. 

. ForK <kP{P),£p,k = {0}. 

• For K > kP {(3), there exists z = z{f3, K) > such that Sp^K = {±z}. 

The next theorem gives the type of each of these points in Sp^x- The type is always 1 except when K = 
Kc^^ {(3); in this case the global minimum point at has type r = 2 if /3 < /3c and type r = 3 if = /3c. 

Theorem 6.1. Let (3 < (3c = log 4 and define Kc {(3) by (6.2). The following conclusions hold. 
(a) For K < kP (/?), z = has type r = 1. 
{h)LetK = kP{P). 

(i) For (3 < Pc< z = has type r = 2. 

(ii) For (3 = (3c, z = {) has type r = 3. 

(c) For K > Kc {(3) and each choice of sign, z = ±z(/3, K) has type r = 1. 

Proof, (a) By (6.2), we have 

G;;,^(0) = 2(3K{l-2(3Kd'pm 
( K \ 
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Therefore, K < kP{P) imphes that Cp xiO) > and thus that z = has type r = 1. 

(b) For K = kP{I3), G^_^(0) = 0. A simple calculation yields 

Therefore, for {3 < (5c, G^''^(0) > and for (3 = Pc, G^p^^^i^) = 0. Computing the sixth derivative yields 

Gg)^(O) = 2{2pKf/9. (6.4) 

As a result, z = has type 2 if /3 < /J^ and has type 3 if P = (3c- 

(c) To prove that the symmetric minimum points ±5 of Gp^K each have type r = 1, we employ the 
results of Lemma 2.3.5 in [18]. This lemma states the existence and uniqueness of nonzero global minimum 
points ±.w = ±w{P, K) of Fp^x{w) = uP' /A(3K — cp{w); iZ; is a global minimum point of Fp^j^ if and 
only if z = w/2(3K is a global minimum point of Gp^K- Lemma 2.3.5 in [18] also states that F'^ j^{w) > 0. 
Since Fp j^{w) > if and only if Gp j^{z) > 0, the symmetry of Gp^K allows us to conclude that for each 
choice of sign ±5 has type r = 1. This completes the proof. ■ 

We next classify by type the points in Sp^K for (3 > (3c and K > 0. According to Theorem 3.7, there 
exists a critical value K^^ with the following properties. 

. ¥or K <K';^\(3),Ep,K = m. 

• For K = i^i^^ {(3), there exists z = z{(3, K) > such that Sp^x = {0, ±z}. 

. For K>K^^\(3),£0,K = {±z}. 
The next theorem shows that the type of each of these points in Sp^x is 1. 
Theorem 6.2. Let (3 > (3c and K > Q. The points in Sp^K all have type r = 1. 

Proof. We first consider when £p^K contains 0, in which case K < K^^\p). Define K2 = l/2Pc'^{0). It 
is proved in Lemma 2.3.15 in [18] that for (3 > (3c^^ have i^c^^ {(3) < K2. Since 

G"p^k{^) = 2(3K{1 - 2(3K4{0)) (6.5) 

it follows that whenever K < K^^\p), 1 > K/K2 and thus Gp ^^{O) > 0. We conclude that the minimum 
point of Gp^K at z = has type r = 1, as claimed. 

For K > Ki^\p), £p^K also contains the symmetric, nonzero minimum points ±z of Gp^K- We prove 
that each of these points has type r = 1, employing the results of Lemma 2.3.11 in [18]. This lemma 
states the existence and uniqueness of nonzero global minimum points ±w = ±w{(3, K) of the function 
Fp^x{w) = /4:(3K — cg{w); it; is a global minimum point of Fp^x if and only if z = w/2(3K is a global 
minimum point of Gp^x- Furthermore, Lemma 2.3.11 states that F^ ^.{w) > 0. Since F'^ j^{w) > if and 
only if G'^ p^{z) > 0, the symmetry of Gp^K allows us to conclude that for each choice of sign ±z has type 
r = 1. This completes the proof. ■ 

We now state the limit theorems for the P„^^^;^-distributions for appropriately scaled partial sums Sn- 
The first. Theorem 6.3, states hmit theorems that are vahd when Gp^K has a unique global minimum point 
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at z = 0. This is the case for (5 < Pc, K < kP{P) [Thm. 3.6(a)] and for (3 > f3c, K < K^\0) [Thm. 
3.7(a)]. The second, Theorem 6.4, states conditioned Umit theorems that are valid when G/3,a- has multiple 
global minimum points. Because Theorems 6.3 and 6.4 are immediate applications of Theorems 2.1 of [9] 
and 2.4 of [10], respectively, we state them here without proof. 

In Theorem 6.3 /q denotes the density of an iV(0, a^) random variable with 

. 2/?frcg(0) 

When the type of the minimum point at is r = 1, cr^ > because in this case Gp ^(0) > and in general 
c^(0) > 0. If / is a noimegative, integrable function on IR, then for r = 1, 2, or 3 we write 

Pn,pMSn/n^-^/^'- G dx} =^ fix) dx 

to mean that as ra ^ 00 the P„,/3,i^-distributions of Sn/n^~^^'^'^ converge weakly to a distribution with 
density proportional to /. When r = 1, / = fQ^„2, and the limit is a central-limit-type theorem with 
scaling n^/^. When r = 2 or 3, the limits involve the nonclassical scaling n^/^ or rc'l^, respectively, and 
the P^ g j^ -distributions of the scaled random variables converge weakly to a distribution having a density 
proportional to exp[— const • z^] or exp[— const • z^]. 

Theorem 6.3. Suppose that Sp^K = {0} and let r he the type of the point z = Q as given in Theorems 6.1 
and 6.2. With the positive quantity defined in (6.6), as n ^ 00 

P / '^^ (zrl \ ^ I forr = l 

"'^'■^lni-V2r ^axj^<^ exp(-Gg2(0)a;2V(2r)!) dx forr = 2orr = 3. 

When r = 2[K = Kf\f3), (3 < Pel G^%(0) is given by (6.3), and when r = 3[K = kP{P), p = /3c], 
Ggi,(0) = 2(2/3i^)V9. 

Our last theorem states conditioned limit theorems that are vaUd when G^^k has multiple minimum 

(2) 

points. This holds in three cases: (1) when (3 < (3c and K > K}- {(3), in which case the minimum points 
are ±z{j3,K) with z{l3,K) > [Thm. 3.6(b)]; (2) when (3 = (3^ and K = K^^\p), in which case the 
minimum points are 0, ±z{l3, K) with z{(3, K) > [Thm. 3.7(b)]; (3) when p = PcWdK > (/?), in 
which case the minimum points are ±z(/9, K) with z{f3, K) > [Thm. 3.7(c)]. In each case in which Gp^K 
has multiple minimum points, Theorems 6. land 6.2 states that all the minimum points have type r = 1. 
Hence, if we denote the minimum points by Zj for j = 1,2 or for j = 1, 2, 3, then for each j we have 
G"p j^{zj) > 0. Defining 

, . 2/?K4(2/3Kz,-) 
" — G" (z ) — ' ^ ^ ^ 

we see that > 0. 

Theorem 6.4. Suppose that £p,K = {^i) ■•■) Zj^} for m = 2 or 3. For each j = 1, . . . ,m we let fg ^2 be 

the density of an N{0, cr|) random variable, where crj is the positive quantity defined in (6.7). Then there 
exists A = A{zj) > such that for any a G (0, A) 

( S — nz- S 1 

Pn,l3,K\ "^^1/2 ^ ~^ ^ ~ '^^\\ =^ fo,aj (x) dx aS U ->■ OO. 

This completes our study of the limits for the Pji^^^^-distributions of appropriately scaled partial sums 
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